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Resonant X-ray diffraction studies on the charge ordering in magnetite 
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Here we show that the low temperature phase of magnetite is associated with an effective, although 
fractional, ordering of the charge. Evidence and a quantitative evaluation of the atomic charges are 
achieved by using resonant x-ray diffraction (RXD) experiments whose results are further analyzed 
with the help of ab initio calculations of the scattering factors involved. By confirming the results 
obtained from X-ray crystallography we have shown that RXD is able to probe quantitatively the 
electronic structure in very complex oxides, whose importance covers a wide domain of applications. 

PACS numbers: 71.30.-|-h 61.10.Nz 78.70.Ck 
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Known in ancient times as lodestone and used to mag- 
netize the mariner's compass 0, magnetite (Fe304) is 
still today the archetype compound of a number of phys- 
ical properties and applications. For instance, it is a 
promising candidate for the development of highly sen- 
sitive magneto-resistive devices in spin electronics. A 
sufficient comprehension of magnetite is mandatory in 
view of its applications and as a reference for similar 
effects in related materials. Particularly, the magnetic 
and magneto-electric properties of magnetite exhibit a 
conspicuous anomaly at Ty — 121K that still defies un- 
derstanding 0. First suggested by Verwey and ex- 
pected in many oxides the low temperature phase 
transition in magnetite has been associated to a charge 
disproportion on the metal atoms sites although a direct 
confirmation has never been evidenced in magnetite up 
to now. Theoretical predictions ^ i3] have recently sup- 
ported Verwey's scheme of the localization of the charge 
(but on a more complex pattern) and nuclear magnetic 
resonance 8] and Mossbauer j£j experiments are compat- 
ible with different oxidation states of the octahedral iron 
sites. However direct confirmation of charge dispropor- 
tionation is still lacking. 

From the structural point of view the determination 
of the atomic positions issued from the metal-insulator 
transition challenges the scientific community ever since 
Verwey's seminal work Progress in the solution of 

the problem runs parallel to the development of new and 
sophisticated experimental techniques as well as to the 
implementation of refined computing codes and appro- 
priate data analysis strategies. Despite all these technical 
advances, magnetite still remains a rather difficult case 
for conventional crystallography: the symmetry lowering 
{Fd3m — > Cc, and Oc x ac x Oc \/2ac x \/2ac x 2ac, 
with ac=8.394 A) generates 8 and 16 non equivalent 
iron sites at tetrahedral and octahedral positions, respec- 
tively, each one with its own atomic charge. However, the 
final structure is not yet perfectly known. The actually 
best refinement has been recently performed by Wright 



and collaborators [10| (in space group Pmca, with lat- 
tice parameters « ac/V^ x Uc/V^cLc x 2ac). Complexity 
is greatly reduced in this structure model (Fig. where 
there are 6 non equivalent iron atoms, two in tetrahedral 
sites (Fet) and four in octahedral sites (Fei, Fe2, Fea and 
Fe4). Only the octahedral irons are supposed to undergo 
a charge ordering (CO). Whereas Fei and Fe2 are at the 
centers of their respective octahedra, Fe^ and Fe4, posses 
more distorted local environments, the iron atoms are 
slightly off center. It can also be noticed that these po- 
sition shifts are responsible for the doubling of the cell 
along the c— axis whereas the Fci and Fe2 ordering has 
no effect on the cell doubled. 

Resonant X-ray diffraction is a technique where both 
the power of site selective diffraction and the power of 
local absorption spectroscopy regarding atomic species 
are combined to the best [llj. Reflections are recorded 
over some tens electron-volts around the absorption edge 
of the element(s) present in the material, where they 
show strong energy and angular dependencies. This 
phenomenon is due to the virtual photon absorption- 
emission associated with the resonant transition of an 
electron from a core level to some intermediate state 
above the Fermi level. By virtue of the dependence on 
the core level state energy and the three dimensional elec- 
tronic structure of the intermediate state, this technique 
is specially suited to study charge orbital 0,0] 

or spin orderings and associated geometrical distor- 
tions 0, 01 . In the case of charge ordering, we exploit 
the fact that atoms with closely related site symmetries 
but with barely different charges exhibit resonances at 
slightly different energies. Here we show that the sen- 
sitivity of this effect allows for quantitative estimations 
of the charge disproportion. Opposite to fluorescence or 
absorption measurements, the power of diffraction relies 
on the capability of detecting differences that are even 
smaller than the inverse lifetime of the core hole level. 
Clearly, not all Bragg reflections are sensitive to charge 
ordering, and the failure in detecting CO signatures in 
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FIG. 1: (Color online) (Left) A sketch of the low tempera- 
ture Fe304 structure in the orthorhombic space group Pmca. 
There are six non equivalent iron atoms, two in the tetrahe- 
dral sites (Fet) and four in the octahedral sites (Fei, Fe2, Fes, 
Fe4) and it is understood that charge ordering occurs at the 
octahedral sites. The actual structure is known to be Cc in a 
approximately %/2fflc x \/2ac x 2ac unit cell, with 16 unequiv- 
alent Fe-sites at octahedral positions. However, no structural 
refinement in this space group has been successful yet. (Right) 
Cell doubling along a (|(oc + he)) and b (|(ac — 6c)) axis of 
the Pmca structure. The color code tags represents Fes and 
Fe4 positions in the settings of the Cc space group. This CO 
model has been tested and greatly improves the quality of the 
refinement. 



previous RXD work on ma gne tite is partly related to an 
incomplete choice of them jlfll E(l | . 

Our RXD experiments were carried out at the U.K. 
CRG beamline (XmaS) in the European Synchrotron Ra- 
diation Facility (Grenoble) operating at photon energies 
between 7060 and 7220 eV. The X-ray beam wavelength 
was selected by a Si(lll) double crystal monochromator 
with an energy resolution of 0.8 eV at the photon en- 
ergies of the experiment. After scattering an analyzer 
crystal was used to filter out background coming from 
the sample (mostly fluorescence). The two samples used 
have been prepared using the vertical floating zone tech- 
nique, with the help of an optical furnace (Crystal System 
Inc. FZ-T-IOOOO-H-VI-P-G). The Fe304 feed rods were 
prepared from high purity powder (Alfa Aesar 99,997%), 
compacted into cylindrical shape under isostatic press 
(up to 1200 Bar). The Verwey transition was charac- 
terized by SQUID measurements and results have been 
given m ref. Vm The as-grown crystals exhibit a sharp 
bulk Verwey transition with Ty = 121 K. The scattering 
experiment was carried out in reflexion mode, off a {110} 



surface for both samples and at 50K. 
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FIG. 2: (Color online) (Top) Energy variation of the intensi- 
ties at the peak for a number of Bragg refiections around the 
Fe K-edge and at T=50K. From top to bottom the reflections 
are: (-4 4 1), (-4 4 2), (-4 4 3), (-4 4 4), (-4 4 5) and (-3 3 3/2). 
Reflections are indexed according to the Miller indices of the 
high temperature cubic phase. Data are completely repro- 
ducible between the two studied samples. The small peaks at 
around 7070 eV in the scans of the (-4 4 1) and (-3 3 1.5) re- 
flections are multiple scattering events that can be shifted in 
energy by a slight change in the azimuth angle. The log scale 
allows to appreciate the quality of the data, mainly above the 
edge. (Bottom) The same reflections corrected from absorp- 
tion. Actual intensities have been normalized so as to fit into 
the graph. 

RXD experiments must be performed on a single crys- 
tal, if a good signal-to-noise ratio is required. Tradi- 
tionally one of the major difficulties of crystallographic 
work on magnetite has been to overcome a very strong 
self absorption, twinning, extinction and the presence 
of multiple-scattering. These issues can be easier over- 
come in RXD studies than in conventional single crystal 
diffraction experiments. Indeed, parameter refinement in 
the former is carried through the analysis of the energy 
dependence of the photons scattered within each of the 
measured Bragg reflections, the absolute scale being a 
unrefined proportionality constant. In order to minimize 
the number of crystallographic domains, a magnetic field 
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of 0.3T has been applied parallel to one of the [001] di- 
rections that uniquely defines the c— direction. Multiple- 
scattering events are eliminated by the measurement of 
the energy spectra at different azimuths and crystal set- 
tings. Raw data have been plotted in figure 13 top. Fi- 
nally, and more important, the absorption correction (fig- 
ure|21 bottom) is carefully performed by using the energy 
dependence of the intensities of several very strong Bragg 
reflections as a measure of absorption. The anomalous 
contribution of these reflections is relatively small and 
can be easily computed. 

To account for the uncertainty of the crystallographic 
structure and the fact that the charge ordering must be 
disentangled from the associated atomic displacements, 
a complete methodology was developed. It needs a very 
important set of experimental data, first principle sim- 
ulations and the use of objective confidence factors for 
comparing experiment and theory |2^. l23j|. Thus, no less 
than 50 reflections, were recorded over an energy range of 
120 eV around the Fe K-edge and at T=50K. This collec- 
tion of peaks is at the variance of the RXD data collected 
in rcfs. 19, 20] and use to dismiss CO in Fe304. 32 reflec- 
tions were kept for the optimization of the parameters, 
12 were reflections indexed in the Cc cell (not used in this 
reflnement) and 6 strong reflections were used for normal- 
ization purposes. The calculations are performed using 
the state of the art code FDMNES jUlil. It calculates 
the flnal states and the resulting structure factors with 
first principal theories. In the present study, the multiple 
scattering theory was used. This code has already been 
efficient in the simulation of different spin, orbital, charge 
or geometrical ordering phenomena ,2^ j27j . 

Trial calculations were performed for the two low tem- 
perature phase structures proposed by lizumi et al. [l^ . 
As expected, calculations in the non centrosymmetric 
space group (Ptoc2i) produced a modest agreement with 
the experimental data (not shown). Calculations in the 
Pmca structure produces a far better agreement and 
in the following we have used the recent results of ref. 
[loj with the charge occupancies on the 3d iron sites as 
the only adjustable parameters. From the local struc- 
tures we conjecture that charge ordering can be charac- 
terized by two different charge disproportions, S12 and 
(534, between the centered Fei (charge state represented 
as Fe+(2-^-''i^)) and Fes (Fe+^^-^+'^i^^) sites and the off- 
centered Fes (Fe+(2-5-'534)) and Fe4 (Fe+^^-^+'^a-*)) sites, 
respectively. Tests on other orderings, including the 
tetrahedrally-iron coordinated or other combinations of 
octahedral iron atoms, have systematically yielded rather 
poor confidence facto rs. 

The most important result is that a significant im- 
provement of the data fitting is achieved when an explicit 
CO is introduced (Fig. ISJ. Reflections such as (-4 4 1) 
are highly sensitive to the magnitude of 612- All of them 
exhibit a double structure around 7126-7131 eV which 
can be taken as a signature of the charge ordering phe- 
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FIG. 3: (Color online) Some of the experimental and cal- 
culated resonant X-ray diffraction peak spectra in Fe304 at 
T=50K. Experimental data have been corrected for absorp- 
tion. Shown are the experimental data (black dots), the cal- 
culated spectra with (red line) and without (blue dots) charge 
ordering. Some reflections such as the (-1 1 0) and the (-4 4 1) 
clearly display the derivative effect (see text) and are therefore 
very sensitive to the CO. Some others, as the (-4 4 2), exhibit 
a line shape that does not arise from the charge ordering. It 
is an extinct Bragg reflection of cubic phase that displays fea- 
tures related to the anisotropy of the Fe-form factors[l^.ll8|. 
Strong reflections, such as the (-4 4 4), were used for normal- 
ization purposes. Half-integer reflections, such as (-4 4 5/2), 
are sensitive to the charge difference between Fes and Fe4. Re- 
flections are indexed according to the high temperature cubic 
phase settings. 



nomenon as seen by RXD. They are a clear example of 
the derivative line shape, characteristic of CO as like- 
wise observed in other compounds 0,0. Indeed, for 
a specific charge pattern and the reflections under con- 
sideration, the partial structure factor associated with 
oppositely charged iron atoms almost cancel. The opti- 
mization of the parameters at the Fei-Fe2 sites gives a 
(5i2 = 0.12±0.025 electrons which corresponds to a shift 
of the absorption edge by 0.9 eV. This result is particu- 
larly robust and remains unchanged whether a large or a 
small data basis is used in the refinement. The effect of 
(^34 is contained in half-integer reflections and even more 
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FIG. 4: (Color online) Iso-value of the confidence factor, the 
metric Di of Philip and Rundgren 22]. The best agreement is 
obtained for 5i2 = 0.12±0.025 electrons and ^34 = 0.10±0.06 
electrons. It corresponds to the 5.38, 5.62 and 5.40, 5.60 
charge occupancies of Fei to Fe4, respectively. 



in reflections indexed in the Cc cell. Interestingly, and 
despite of the lack of knowledge about the actual lattice 
distortion yielding a Cc unit cell, the charge can be re- 
fined in this space group (see Fig. ^ right, for the charge 
distribution compatible with the Cc symmetry) by using 
reflections indexed for the Pmca symmetry. The value 
of (534 is associated with a large uncertainty: 0.10±0.06 
electrons, most probably reflecting the fact that each of 
the Fes and Fe4 positions gives rise to four noncquiv- 
alent sites, being not considered in this work. To ob- 
tain a precise estimate of would require a concomi- 
tant knowledge of the atomic displacements producing 
the Cc structure and the use of reflections indexed in the 
\/2ac X \/2ac x 2ac unit cell. Figure 0] represents the con- 
fidence factor plotted in the form of a contour plot. All 
the refined charges are in very good agreement with the 
theoretical predictions using LDA-I-U calculations 0, Q • 
The contours run nearly parallel to the 5^^ axis mean- 
ing that there is little correlation between both charge 
orderings. 

In this work we have solved the long standing contro- 
versy of whether or not actual charge ordering develops at 
the octahedral Fe-sites of Fe304 as a result of the metal- 
to-insulator (Verwey) phase transition. Several conclu- 
sions can be drawn out. First, we have shown that signif- 
icant, non- integer charge localization occurs at the metal 
sites in the low temperature phase of magnetite. This re- 
sult is robust and independent of the uncertainties of the 
actual crystallographic structure. In view of the short 
interaction time of the photons, « 10~^^ s., we conclude 
that the ordering is purely static. Second, the results of 
our experiment and additional calculations considering 
the magnitudes of the charge are in good overall agree- 
ment with those resulting from the bond valence sum 
method (BVS). Third, CO at Fei and Fe2 sites gives rise 
to a [001] c charge modulation whereas CO at Fe3 and 
Fe4 originate the [00^]c ordering, and both are relatively 



uncoupled. Fourth, the magnitude of the charge is in 
agreement with the several band structure calculations 
[a, Q which implies that the phase transition, at least 
the ordered part, can be explained under electronic con- 
siderations alone. 

The authors want to thank V. Dmitrienko and S. 
Fratini for helpful discussions and A. Kirfel for a criti- 
cal reading of the manuscript. 

Note Added.- While the mansucri pt w as with referees 
we became aware of a RXD paper |22| where evidence 
for the CO has been found, thus supporting our results. 
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